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Abstract 

We have studied the electronic structure of the Ni triangular lattice in NiGa2S4 using photoe- 
mission spectroscopy and subsequent model calculations. The cluster-model analysis of the Ni 2p 
core- level spectrum shows that the S 3p to Ni 3d charge-transfer energy is ~ -1 eV and the ground 
state is dominated by the d 9 L configuration (L is a S 3p hole). Cell perturbation analysis for the 
MS2 triangular lattice indicates that the strong S 3p hole character of the ground state provides 
the enhanced superexchange interaction between the third nearest neighbor sites. 

PACS numbers: 75.30.Et,79.60.-i,75.50.-y 
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Unusual magnetic properties of geometrically frustrated spin systems have been attract- 
ing broad interest of many theorists and experimentalists in the field of condensed matter 
physics In frustrated spin systems with orbital degeneracy, specific orbital patterns can 
lift the magnetic frustration and some long range magnetic orderings can be realized in the 
ground state . In frustrated spin systems without orbital degeneracy, spin-lattice 

coupling tends to provide cooperative distortion lifting the magnetic frustration 0, Is, or 
small randomness tends to induce spin freezing at low temperature [k]]. Therefore, although 
spin disordered states including the resonating-valence bond state are expected in the ge- 
ometrically frustrated spin systems [III, |l2[], few systems show spin disordered states as 
ground states 13|. Among the various magnetic compounds with geometrically frustrated 
lattice, the newly-discovered NiGa2S4 has the Ni 2+ (S = l) triangular lattice layer without 
orbital degeneracy (see Fig. [1]) and is found to have a frozen spin-disordered ground state 
using neutron diffraction and nuclear quadrupole resonance experiments [lj], Ijj]. The neu- 
tron result also indicates that the spin-spin correlation between the third nearest neighbors 
is much stronger than that between the first and second nearest neighbors, indicating that 
the conventional triangular lattice model with the nearest neighbor superexchange coupling 
is not enough to describe NiGa 2 S4 jl6l |. 

In order to understand the origin of the unusual magnetic properties of NiGa2S4, it is 
highly important to reveal its underlying electronic structure. In this Letter, we report 
photoemission study of NiGa2S4 single crystals. The photoemission results and subsequent 
model calculations indicate that the ground state has the d 9 L character (L is a S 3p hole) and 
that the strong S 3p hole character of the ground state provides the enhanced superexchange 
interaction between the third nearest neighbor sites. NiGa 2 S 4 is a unique spin-disordered 
system in that the negative charge-transfer energy allows relatively long superexchange 
pathways. 

Single crystals of NiGa2S4 are grown by chemical vapor transport as described in ref. 
17| . The x-ray photoemission spectroscopy measurements were performed using a JPS 9200 
spectrometer equipped with a monochromatized Al Ka x-ray source {hv = 1486.6 eV). The 
total energy resolution was ~ 0.6 eV. The single crystals were cleaved in situ in order to 
obtain clean surfaces. All photoemission data were collected at room temperature within 
48 hours after cleaving. As expected from the layered structure of NiGa2S4, the obtained 
surfaces were extremely clean and stable: the photoemission spectra were consistent with 
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FIG. 1: Schematic drawing for the Ni 2+ (5=1) triangular lattice layer and the electronic config- 
urations of the NiS6 cluster model used to analyze the photoemission spectra. The N1S2 layer is 
constructed from the NiSg octahedra in which one Ni ion is located at the center and six S ions at 
the verteces. In the NiSg cluster, the S 3p molecular orbitals with e 9 and t2 g symmetry hybridize 
with the Ni 3d orbitals with e g and t2 g symmetry. L denotes a hole in the S 3p molecular orbitals. 
The charge-transfer energy A is given by the excitation energy from d s to d?L. 





1 1 1 1 1 1 

NiGa 2 S 4 






S 2p XPS 








NiS, | W^/ GaS 


Intern 




|2p 3/2 , ft 
jg> / Ga 3s^^to,«^^ 




1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 




168 166 


164 162 160 



Binding Energy (eV) 

FIG. 2: S 2p and Ga 3s core-level XPS of NiGa2S4 (open circles). The spectrum is fitted to five 
Gaussians: the Ga 3s component (thin solid curve), the 2p 3 ^ 2 and 2pi/ 2 components of the GaS 
layer (dashed curve), and those of the NiS2 layer (dotted curve). The fitted result is shown by the 
thick solid curve. 

the stoichiometric surface and did not show any change during the experiment. 

The S 2p core-level photoemission spectrum of NiGa2S4 is shown in Fig. [2j The S 2p 
spectrum can be decomposed into four components: the 2p 3 / 2 and 2pu 2 components of the 
GaS layer and those of the N1S2 layer. The S 2p peaks of the MS2 layer are higher in 




FIG. 3: Ni 2p core-level XPS of NiGa2S4 (open circles). The calculated line spectrum is broadened 
(solid curve) and is compared with the experimental result. In the lower panel, the line spectrum 
is decomposed into the cd s , cdrL, and cd 10 L 2 components. 



binding energy than those of the GaS layer. The energy difference of ~ 1 eV indicates that 
the S ions in the MS2 layer have less valence (S 3p) electrons than those in the GaS layer. 
This observation already suggests that the amount of the S 3p holes in the N1S2 layer is 
substantial as revealed by the Ni 2p spectrum in the next paragraph. 

The Ni 2p core- level spectrum of NiGa2S4 is shown in Fig. [3j The Ni 2p 3 / 2 spectrum 
consists of three structures: the main peak at ~ 854 eV and the satellite structures at ~ 
859 eV and 864 eV. In order to extract the electronic structure parameters such as the S 
3p to Ni 3d charge-transfer energy A, the Coulomb interaction between the Ni 3d electrons 
U, and the transfer integrals between the S 3p and Ni 3d orbitals (pda), we have performed 
configuration-interaction calculations using an octahedral NiSg cluster model (see Fig. CD). 
Since the NiS6 octahedra form the MS2 layer sharing their edges, the interaction between 
the neighboring NiSg octahedra via the S 3p a orbitals pointing to the Ni sites is rather weak 



compared to that in the corner sharing case [181 ] . Therefore, the cluster-model analysis is 
expected to give a good description of the Ni 2p core-level spectrum. In the present cluster 
model, the Coulomb interaction between the Ni 3c? electrons are given by the Slater integrals 
F°(3d, 3d), F 2 (3d,3d), and F A (3d, 3d). The average Ni 3c?-Ni 3d Coulomb interaction U is 
expressed by F°(3d, 3d) and is an adjustable parameter. F 2 (3d, 3d) and F A (3d, 3d) are fixed 
to 80% of the atomic Hartree-Fock values [19( . The Coulomb interaction between the Ni 2p 
core hole and the Ni 3d electron is expressed by the Slater integrals F°(2p, 3d), F 2 {2p, 3d), 
and G l {2p, 3d). The avera ge N i 2p-Ni 3c? Coulomb interaction Q is expressed by F°(2p, 3d) 
and is fixed to ^/0.8_|2o|, \2^. F 2 (2p,3d) and G\2p,3d) are fixed to 80% of the atomic 
Hartree-Fock values [191 ]. 

In Fig. [31 the calculated line spectrum is broadened and compared with the experimental 
result. The three-peak structure of the Ni 2p^/2 spectrum is well reproduced by the calcu- 
lation using A=-1.0 eV, {7=5.0 eV, and (pda)=-1.0 eV. The ground state is given by the 
linear combination of d 8 , d 9 L, and d 10 L 2 configurations: 

^ g = a\d 8 )+(3\d 9 L}+^\d w L 2 }, (1) 

where L denotes a hole in the S 3p orbitals. The final states are given by the linear combina- 
tions of cd s , cd 9 L, and cd 10 L 2 configurations, where c denotes a Ni 2p core hole. In the final 
states, the Coulomb interaction between the Ni 2p core hole and Ni 3c? electrons stabilizes 
cd l0 L 2 compared with cc? 9 L. The present analysis gives a 2 = 0.25, f3 2 = 0.60, and 7 2 = 0.15, 
and the ground state is dominated by the d 9 L configuration. Negative charge-transfer en- 
ergy compounds have been found in some transition-metal oxides with high valence such 



as Cu 3+ and Fe 4+ 22]. The present system is the first example of transition-metal sulfide 
with negative charge-transfer energy. The negative charge-transfer energy indicates that the 
ground state has strong S 3p hole character. Although the symmetry of the ground state is 
3 A 2g , the magnitude of the local spin at the Ni site is reduced to 0.55 in comparison with 
the ionic value S — 1 by the strong Ni 3c?-S 3p hybridization. In this sense, the 3 A 2g state 
with the d 9 L character is a triplet analogue of the Zhang-Rice singlet in high-Tc cuprates. 

In Fig. HI the valence-band spectrum is plotted with the line spectrum obtained by 
the NiSe cluster model calculation with the parameters obtained from the Ni 2p spectrum. 
Although the contribution from the GaS layer prevents us from fitting the valence-band 
spectrum, the overall agreement between the line spectrum and the valence-band spectrum 
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FIG. 4: Valence-band XPS of NiGa2S4 (open circles). The line spectrum is obtained by the NiSq 
cluster model calculation. 



support the cluster-model analysis of the Ni 2p spectrum. The first ionization state has the 
symmetry of 2 E g and the second ionization state is 4 T l9 . This situation is similar to those of 



26] 



NiS and N1S2 which have small but positive charge-transfer energy [24J, [25|]. The excitation 
from the 3 A 2g ground state to the 2 E g first ionization state is obtained by removal of an e g 
electron. The 2 E g state is expected to have some energy-momentum dispersion due to the 
interaction between the NiSq clusters and to form the e g band. 

In the next step, let us examine superexchange pathways in NiGa2S4 based on the elec- 
tronic structure parameters obtained from the analysis of the Ni 2p spectrum. Since the t 2g 
orbitals are fully occupied and the neighboring Ni-S bonds are approximately orthogonal to 
each other in the two NiS6 clusters sharing their edges, the cell perturbation calculation 
becomes rather simple and the superexchange interaction between the two clusters is given 
by the second order perturbation of S 3p-3p transfer term. In the NiS 6 cluster, the exact 
ground state is given by ^ g ( 3 A 2g ) = a\d 8 ( 3 A 2g )} + (3\d 9 L( 3 A 2g )} + -f\d 10 L 2 ( 3 A 2g )}. Here, the 
NiSe cluster is assumed to be a cubic octahedron in the calculation although it has a slight 
trigonal distortion in NiGa 2 S 4 . Since the S 3p hole denoted as L is dominated by the Sub- 
orbitals pointing to the Ni sites, it is reasonable to approximate that 

\d 9 L( 3 A 2g )) = —(\d x 2_ y 2L 3z 2_ r 2 \ + \L x 2_ y 2d 3z 2_ r 2\). (2) 

Here, d x 2_ y 2 and d 3z 2- r 2 denote a hole in Ni 3d x 2 — y 2 and 3z 2 — r 2 orbitals, respectively. 




FIG. 5: (Color online) (a) Top view of the nearest neighbor clusters in the MS2 layer, (b) Top and 
side views for superexchange pathway bewteen the second nearest neighbor clusters, which is given 
by the S 3p-3p transfer of — (ppir), as indicated by the thick solid line, (c) Top and side views for 
two superexchange pathways bewteen the third nearest neighbor clusters, which are given by the 
S 3p-3p transfer of [— (ppu) + (jrpir)]/y/2, as indicated by the thick solid and broken lines. 

The x-, y-, and z-axes are along the Ni-S bonds as shown in Fig. [5](a). The x 2 — y 2 type S 3p 
molecular orbital is given by L x 2_ y 2 = \{p x — p 2 +p x — p y ), in which pi, p 2 , p x , and p y are p x 
and p y orbitals pointing to the Ni site [see Fig. [3(a)]. The d 3z 2_ r 2 type S 3p molecular orbital 
is given by L 3z 2__ r 2 = -^(p 3 z +p 6 z ) — -^(pl +p 2 y +P X +P y ), in which p 3 z and p\ are p z orbitals 
pointing to the Ni site. As a consequence, the overlap integral between the neighboring 
d 9 L states with 3 A2 9 symmetry can be neglected. The superexchange pathways between 
the neighboring <i 9 L( 3 A 2s ) states are given by the transfer terms between the two S 3p 



molecular orbitals with L x 2_ y 2 and L 3z 2_ r 2. The transfer integrals between the neighboring 
L x 2_ y 2 orbitals is \[{ppcr) + (pp7r)] and that between the neighboring L 3z 2_ r 2 orbitals is 
12 \{PP a ) +9(PP 7I ")]5 where (ppa) and (ppn) are the S 3p-3p transfer integrals. Therefore, the 
antiferromagntic superexchange interaction between the neighboring sites is given by 

where U p is the Coulomb interaction between the S 3p holes at the same S sites and is 
expected to be ~ 1 eV 0. Since (3 2 = 0.6, (ppa) = 0.6 eV, and {ppn) = -0.15 eV, J^ F ~ -5 
meV. On the other hand, the two S 3p holes at the shared sulfur sites have Hund coupling 
due to the on-site exchange integral J p at the S site. The ferromagnetic interaction due to 
the Hund coupling is given by J[ = j^J p ~ 4 meV, where J p is the Hund coupling between 
the S 3p holes at the same S sites and is assumed to be ~ 0.2 eV. As a result, the total 
magnetic interaction between the neighboring sites J\ = J^ F + Jf is estimate to be ~ -1 
meV. The transfer integral between the second neighbor L x 2_ y 2 orbitals is given by — \{ppTi) 
and that between the second neighbor L 3z 2_ r 2 orbitals can be neglected [see Fig. [5] (b)]. 
The superexchange interaction J 2 due to the transfer term between the second neighbor d 9 L 
states becomes 

which is ~ -1 meV and as small as J\. 

In contrast, the transfer integral between the third neighbor L x 2_ y 2 orbitals is \[{ppo) — 
(ppn)] and that between the third neighboring L 3z 2_ r 2 orbitals is {pP a ) + {PP 17 )] as 

shown in Fig. [5](c). The superexchange interaction J3 between the third neighbor sites due 
to the transfer term between the d 9 L states is given by 

which is estimated to be ~ -14 meV. The superexchange interaction between the forth 
neighbor sites is beyond the second order perturbation on the p-p transfer integrals {ppo) 
and {ppTc), and is expected to be negligibly small. The present cell perturbation analysis 



predicts that the third neigh 



Dor superexchange interaction is dominant consistent with the 



neutron diffraction result 14j 



In conclusion, the electronic structure of NiGa2S4 with the N1S2 triangular lattice is inves- 
tigated using photoemisssion experiment and model calculations. The Ni 2p photoemission 



data and the cluster model calculation show that the ground state has the d 9 L character (L 
is a S 3p hole) which is a triplet analogue of the Zhang-Rice singlet state. The strong S 3p 
hole character of the ground state provides the enhanced superexchange interaction between 
the third nearest neighbor sites. NiGa2S4 is a unique spin-disordered system in that the 
negative charge-transfer energy allows unexpectedly long superexchange pathways. 
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